A COMPARATIVE STUDY OF TWO MODELS OF QGP-FIREBALL 

FORMATION. 
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A Comparative study of the strengths and weakness of the models of 
fireball formation namely the statistical model of Ramanathan et.al (Physical 
Review C 70, 027903, 2004) and the approximation schemes of Kapusta et. al 
(Physical Review D 46, 1379, 1992) and its subsequent improved variants is 
made. The way to complement the various approximation schemes, in order to 
enhance their utility in the phenomenological analysis of QGP data that are 
expected from ongoing URHIC experiments, is suggested. The calculations, 
using the former model, demonstrate a striking QCD behaviour of the surface 
tension of the QGP droplet resulting in its increase with temperature, which 
is due to the confining nature of QCD forces at the surface and the interface 
surface tension varies as the cube of the transition temperature which is in 
conformity with the results of Lattice QCD simulations. 



The formation of QGP droplet (fireball) is one 
of the most exciting possibility in the ultra rel- 
ativists heavy ion collision (URHIC) [1]. The 
physics of such an event is very complicated and 
to extract meaningful results from a rigorous use 
of QCD appropriate to this physical system is 
almost intractable though heroic efforts at lat- 
tice estimation of the problem has been going 
on for quite some time [2]. One way out is 
to replicate the approximation schemes which 
have served as theoretical tools in understand- 
ing equally complicated atomic and nuclear sys- 
tems in atomic and nuclear physics in the context 
of QGP droplet formation. This approach lays 
no claim to rigour or ab-initio "understanding" 
of the phenomenon but lays the framework on 
which more rigorous structures may be built de- 
pending on the phenomenological success of the 
model as and when testable data emerges from 
ongoing experiments. 

In this paper we discuss the strength and weak- 
ness of two approximation schemes which seem 
promising in their usefulness in the above con- 



text and also explore their parametric interre- 
lations which may be useful in phenomenologi- 
cal applications. The two approximation schemes 
under consideration are the QCD oriented QGP- 
droplet formation model of Csernai and Kapusta 
and its later refinements [3] and the simple sta- 
tistical model of Ramanathan et. al [4] which 
again is essentially QCD oriented in that use is 
made of an effective QCD potential. In both 
these schemes it is possible to estimate the fire- 
ball radius, the transition temperatures for their 
formation, nucleation rate etc. in quantitative 
terms for comparision with raw data as and when 
it is available from ongoing URHIC experiments. 

As we shall see the two approaches can com- 
plement each other in the event of analysing 
the data on fireball formation especially with re- 
gard to nucleation of QGP droplets in a hadronic 
medium. 

Central to this comparision is the rate I to 
nucleate droplets of QGP in a hadronic gas per 
unit time per unit volume is given by [5] 
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Where Iq is the nucleation rate at vanishing 
change in free energy AF of the system due to 
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the formation of a single critical size droplet of 
plasma. In the Kapusta et. al [3] formulation it 
is not possible to estimate the whole prefactor in 
terms of other possible measurable parameters of 
the droplet. 

The nucleation process is driven by statistical 
fluctuations being determined by the critical free 
energy difference between two phases. The Ka- 
pusta et. al model [3] uses the liquid drop model 
expansion for this , as given by 
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The first term represents the volume contri- 
bution, the second term is the surface contribu- 
tion where a is the surface tension, and the last 
term is the so called shape contribution . The 
shape contribution is an entropy term on account 
of fluctuations in droplet shape which we may 
ignore in the lowest order approximation. The 
critical radius R c can be obtained by minimising 
(2) with respect to the droplet radius R , which 
in the Linde approximation [6] is, 
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fig. l: The free energy difference AF(R) between 
a hadronic phase with and without a quark- 
gluon plasma droplet (PRC 51, 901; 1995) 
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There is no way to estimate the value of the 
crucial input a within the Kapusta et. al ap- 
proach and the related approaches, though a 
value of 50 MeV/fm 2 is assumed. The criti- 
cal free energy in this model varies as in fig.l [3]. 

It is easy to observe that this is a typical first 
order phase transition behaviour indicated by 
the Ramanathan et.al approach, thereby allow- 
ing direct comparision and the opportunity to 
use the calculational advantages of each of the 
approaches to bolster the usability of the approx- 
imation schemes as a phenomenological tool in 
the analysis of the data as we shall see in the 
following. 

In the approximation schemes of the Ra- 
manathan et. al [4], the relativistic density of 
states for the quarks and gluons is constructed 



adapting the procedures of the Thomas-Fermi 
construction of the electronic density of states 
for complex atoms and the Bethe density of 
states [6] for nucleons in the complex nuclei as 
templates. The expression for the density of 
states for the quarks and gluons (q, g)in this 
model is 
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(4) 

where k is the relativistic four-momentum of 
the quarks and gluons, v is the volume of the 
fireball taken to be a constant in the first ap- 
proximation andV con / is a suitable confining po- 
tential relevant to the current quarks and gluons 
in the QGP [4] given by 



[W( fc )U = (V2*)7 M g 2 (k)T 2 -m 2 /2k (5) 
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where mo is the mass of the quark which we 
take as zero for the up and down quarks and 
150MeV for the strange quarks. The g{k) is the 
QCD running coupling constant given by 

g 2 (k) = (4/3)(127r/27)(l/Jn(l + k 2 /A 2 )) (6) 

where A is the QCD scale taken to be 
150 MeV. 

The model has a natural low energy cut-off at 



with 



k mm = ( 7<?)9 iV 1/3 T 2 A 2 /2) 1 / 4 



N = (4/3)(12tt/27). 



(7) 



The freee energy of the respective case % (quarks, 
gluons, interface etc.) for Fermions and Bosons 
(upper sign or lower sign) can be computed using 
the following expression 
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F t = T Tg z J dk Pi {k) ln(l ± e ^V^)/ T ) (8) 

with the surface free-energy given by a modified 
Weyl [7] expression: 

-^surface ^ 

where the hydrodynamical flow parameter for 
the surface is : 

7 = V2x ^{l/ lg y + {l/ lq )\ (10) 

For the pion which for simplicity represents the 
hadronic medium in which the fireball resides, 
the free energy is : 

F n = (3T/2n 2 )(-v) ^ k 2 dk\n(l-e-^ mi+k2/T ) 
Jo 

(11) 

With these ingredients we can compute the 
free-energy change with respect to both the 
droplet radius and temperature to get a phys- 
ical picture of the fireball formation, the nucle- 
ation rate governing the droplet formation, the 



fig. 2: F total at 7 9 = 67,, 7, = 1/6 for various 
temperatures. 



nature of the phase transition etc. This can be 
done over a whole range of flow-parameter val- 
ues [4] , We exhibit only the two most promising 
scenarios in fig. 2 and fig. 3. 

As could be seen the general nature of the 
curves given by fig. 1 and fig.2 and 3, though us- 
ing two totally different approches leads to the 
conclusion that both predict a first order phase 
transition for the fireball production process. As 
a derivative computation we can also compute 
the nucleation rate leading to the droplet forma- 
tion, which for the two scenarios is illustrated by 
figs. 4 and 5. 

As could easily be observed from figs 4 and 5, 
for both the sets of parameter values, the droplet 
formation rate grows exponentially in the vicin- 
ity of transition temperature range of 150 to 
170 Mev as expected from lattice simulations 
and other model calculations. The critical free- 
energy and radius at the maximum of the curves 
in figs. 2 and 3 allow us to compute the sur- 
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1 : T = 140 MeV 
2:T= 150 MeV 
3: T= 170 MeV 
4: T = 190 MeV 
5: T = 210 MeV 
6: T = 230 MeV 
7 : T = "250 MeV 





FIG. 3: F total at 7 9 
temperatures. 
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fig. 4: The fireball formation rate with tempera- 
tures at 7 3 = 67^, 7g = 1/6 . 



face tension used in the Kapusta et.al and the 
derivative models [3]. 



T c 


AF C 


R c 


a 


a 


(MeV) 


(MeV) 


(fm) 


(MeV/fm 2 ) 


150 


332.203 


3.475 


6.568 


0.078 


160 


382.359 


3.385 


7.966 


0.078 


170 


433.037 


3.285 


9.580 


0.078 


190 


532.219 


3.085 


13.35 


0.078 


210 


623.349 


2.875 


18.004 


0.078 


230 


702.041 


2.655 


23.776 


0.078 


250 


766.041 


2.455 


30.343 


0.078 



Table- 1 for Surface Tension of QGP droplet at 

lg = 879,79 = V 6 - 



T 


AF C 


R c 


a 


a 


(MeV) 


(MeV) 


(fm) 


(MeV/fm 2 ) 




150 


943.595 


5.835 


6.616 


0.078 


160 


1197 


5.965 


8.031 


0.078 


170 


1494 


6.085 


9.633 


0.078 


190 


2216 


6.275 


13.435 


0.078 


210 


3088 


6.375 


18.140 


0.078 


230 


4059 


6.375 


23.844 


0.078 


250 


5052 


6.275 


30.630 


0.078 



Table-2 for Surface Tension of QGP droplet at 

The tables 1 and 2 list the surface tension com- 
puted in the Ramanathan et.al model which can 
be used in the dynamical models[3], thus reduc- 
ing the arbitrariness in the latter models to this 
extent, thus enabling us to use the parameter 
extracted from a static situation to make per- 
dictions about the dynamical growth process of 
fireball formation. 
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fig. 5: The fireball formation rate with tempera- 
tures at 7 g = 87, , 7 g = 1/6 . 

In both tables 1 and 2 the surface tension is 
seen to increase with the temperature of the fire- 
ball, which is a beautiful demonstration of a QCD 
effect. As the temperature of the QGP droplet 
increases the shear forces on the fireball sur- 
face will also increase tending to tear the surface 
quarks apart, consequently, bringing into play 
the confining property of the QCD forces mani- 
festing itself in increased surface tension, which 
is exactly what the calculations show. Another 
striking feature of the result is the independence 
of the QGP droplet surface tension a variations 
in the flow parameters of the model and it varies 
with only the temperature, in the lowest order 
approximation we have employed. 

The constancy of the ratio indicates a cubic 
crtitical temperature dependence of the surface 
tension of the interfacial separation between the 
two phases. This is in striking conformity with 



the results of Lattice QCD simulation [8] . 
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